clc;
clear;
x=[1 2 3 4 5 6 7 8 9 10 ];
y=[20587 12489 5841 3184 1050 244 73 16 12 1 ];
plot(x,y,'r');